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Abstract 

Due to the current trend towards smaller cells, an increasing number of users of cellular networks reside at 
the edge between two cells; these users typically receive poor service as a result of the relatively weak signal and 
strong interference. Coordinated Multi-Point (CoMP) with Joint Transmission (JT) is a cellular networking technique 
allowing multiple Base Stations (BSs) to jointly transmit to a single user. This improves the users’ reception quality 
and facilitates better service to cell-edge users. We consider a CoMP-enabled network, comprised of multiple BSs 
interconnected via a backhaul network. We formulate the OFDMA Joint Scheduling (OJS) problem of determining 
a subframe schedule and deciding if and how to use JT in order to maximize some utility function. We show that 
the OJS problem is NP-hard. We develop optimal and approximation algorithms for specific and general topologies, 
respectively. We consider a time dimension and study a queueing model with packet arrivals in which the service rates 
for each subframe are obtained by solving the OJS problem. We prove that when the problem is formulated with a 
specific utility function and solved optimally in each subframe, the resulting scheduling policy is throughput-optimal. 

Via extensive simulations we show that the bulk of the gains from CoMP with JT can be achieved with low capacity 
backhaul. Moreover, our algorithms distribute the network resources evenly, increasing the inter-cell users’ throughput 
at only a slight cost to the intra-cell users. This is the first step towards a rigorous, network-level understanding of 
the impact of cross-layer scheduling algorithms on CoMP networks. 

Keywords: Coordinated Multi-Point (CoMP), Joint Transmission, Cellular Networks, Approximation Algorithms, 
Scheduling, Queueing Networks. 


I. Introduction 

Cellular networks face an ever-increasing bandwidth demand, driven by the advent of sophisticated mobile devices 
and new applications. Satisfying this demand calls for improvements in the spectral utilization and reductions 
in inter-cell interference. The latter is becoming more relevant as the number of inter-cell users increases with 
ever-decreasing cell sizes. Such users are often unable to receive any transmission due to the high interference. 
Interference reduction can be efficiently accomplished through multi-cell coordination, known as Coordinated Multi- 
Point (CoMP) or Network-MIMO. One particularly promising category is CoMP with Joint Transmission (JT), where 
multiple Base Stations (BSs) jointly transmit to a single user, using the same time-frequency slots. This technique is 
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incorporated in the LTE-Advanced standard [ 11. Recently, CoMP with JT was shown to obtain substantial throughput 
gains in both indoor and outdoor testbeds [20 iQ 

As a result of the implementation of CoMP with JT in the LTE-Advanced standard, algorithm design and 
performance evaluation for these systems have recently received increased attention in the research literature (see, 
e.g„ (T^, (22), (37)). However, most existing work is concerned with developing heuristics designed for saturation 
conditions. In contrast, we consider a cellular network where new packets are generated over time, and construct a 
rigorous framework to develop scheduling algorithms for CoMP with JT that maximize throughput. This is achieved 
via a cross-layer approach, consisting of PHY (considering SINR-based transmission probabilities), MAC (deciding 
on a transmission schedule), and network layer (forwarding traffic over the backhaul). 

We consider a cellular network comprised of multiple BSs interconnected via backhaul links.Users are assigned 
a serving and a secondary BS, and packets destined for a user can be transmitted either by the serving BS only, or 
jointly by the serving and secondary BSs. The latter provides better signal-to-interference-plus-noise ratio (SINR), 
but requires a subframe in both BSs, as well as forwarding the packet from the serving to the secondary BS prior 
to the transmission. A scheduling algorithm for CoMP with JT therefore needs to balance the performance benefits 
of transmitting packets using JT with the additional resources required for doing so. 

We first focus on a single subframe, and study the OFDMA Joint Scheduling (OJS) problem of determining a 
schedule to maximize some utility function, given a set of packets for each user. Such a schedule determines which 
packets should be forwarded over the backhaul and which packets should be transmitted wirelessly, either using 
JT or by the serving BS only. We show that the OJS problem is NP-hard and describe a framework for solving it 
efficiently by decomposing it into problems related to knapsack and coloring. This allows us to develop an efficient 
algorithm for solving the OJS problem in bipartite backhaul network graphs. While backhaul network graphs are 
not necessarily bipartite, this result enables us to develop approximation algorithms for general backhaul graphs. 

We then consider the network evolution over multiple subframes. We define a queueing model where the users 
are fixed, and packets for the various users are generated over time. Departures are determined by the schedule 
obtained from solving the OJS problem in each subframe. We characterize the capacity region (i.e., the packet 
arrival rates that can be sustained). Moreover, we demonstrate that when the OJS problem is formulated with a 
specific queue-length based utility function and solved optimally in each subframe, we obtain a MaxWeight-like 
scheduling policy (e.g., (30)), which we show to be throughput optimal. This is surprising, given the difference 
between OJS and the typical matching-type problems where MaxWeight scheduling performs well. Based on the 
queueing model, we present extensive simulation results to evaluate the performance of the proposed scheduling 
algorithms, as well as the benefits of JT. In particular, we consider different network topologies with an SINR- 
based channel model. We show that the bulk of the gains from CoMP with JT can be achieved with low capacity 
backhaul links. This result is highly relevant as the deployments of advanced cellular wireless technologies have 
a strong impact on mobile backhaul operational expenditure (OPEX), which amount to 20-40% of total mobile 

'There are two flavors of CoMP with JT: coherent j^O] and non-coherent. We consider coherent JT but remark that all results can be directly 
extended to non-coherent JT. 
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operator’s OPEX due to their reliance on Tl/El backhaul copper lines 0- A promising alternative is wireless 
backhaul (e.g., satellite, microwave), which is becoming a viable technology for geographically challenging regions 
and 5G networks. However, since such technology has limited capacity (due to, e.g., limited wireless spectrum and 
poor wireless channel conditions), our result are relevant to both present and future networks. We show that our 
algorithms distribute the network resources more evenly as the backhaul capacity increases. In fact, they increase 
the inter-cell users’ throughput at only a slight cost to intra-cell users. 

The main contributions of this paper are two-fold: (i) we define a rigorous model for CoMP with JT and develop 
novel scheduling algorithms with throughput guarantees for networks with queueing dynamics; (ii) via extensive 
simulations, we observe the benefits of JT and the tradeoffs related to its implementation. 

The rest of the paper is organized as follows. In Section [II] we discuss related work. In Section m we present 
the model. In Section |IV[ we introduce the OJS problem and show that it is NP-hard. We develop approximation 
algorithms for OJS in Section [V] In Section VI we develop and present results for a queueing model, which we 


study through extensive simulation experiments in Section VII Section VIII provides conclusions and directions 
for future research. All proofs can be found in Appendix A. 


II. Related Work 

Previous work on scheduling for CoMP with JT has focused exclusively on analyzing the performance of 
heuristics, and has been limited mostly to networks that are saturated (i.e., have infinite backlog). The proposed 
heuristics are then evaluated via simulations or in testbeds (see, e.g., [20J and references therein). For example, 0 
(22), ( 33 ), 0 , 0 present heuristics for throughput maximization, assuming a perfect backhaul (infinite capacity 
and no delay), while |37| does the same for the finite backhaul case. In [15(, the authors devise a heuristic for 
networks with perfect backhaul, and aim to minimize the backhaul traffic under certain SINR constraints. To our 
knowledge, ©< 0 are the only studies that consider unsaturated networks, where traffic is generated over time 
rather than assuming an infinite supply of available packets. Both of these propose heuristics, assuming a perfect 
backhaul. The main contrast between our work and (6), 0, (22), 0, 0, 0, 0 is that we derive the first 
scheduling policies with performance guarantees. This is done for unsaturated networks, assuming finite backhaul 
capacity and positive delay over the backhaul. 

Models similar to the one considered in this paper have been investigated in the context of single-cell cellular 
transmissions. Packet-level scheduling algorithms for cellular networks are developed in (8), ©• In (3), approxi¬ 
mation algorithms that provide queue stability are analyzed for a single BS. 

Closely-related technologies to CoMP with JT are network-MIMO, multi-cell MIMO, and multi-user MIMO 
(MU-MIMO) @, 0 , 0 , 0 . While theoretical studies (e.g., 0, 0 , 0 ) show that under certain conditions 
such technologies can completely cancel inter-cell interference, achieving these gains in practical scenarios is still 
challenging 0, 0, 0 (e.g., due to the high signal processing complexity). The Study of scheduling schemes 
for these technologies is subject to further research. 
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(a) wireless network 



(b) backhaul network and joint transmission graph 


Fig. 1. A cellular network comprised of basestations (black circles) and users (grey circles). 


III. Network Model 

We consider an OFDMA cellular network comprised of a set of BSs B = {1,..., B} and a fixed set of stationary 


users A f = {l,...,iV}, see Fig. 1(a) The BSs are connected by a backhaul network represented by a graph 
Gj = (B,C), where C is a set of backhaul links with \C\ = C. We refer to Gj as the Joint Transmission Graph , 
as only neighboring BSs in Gj can joint-transmit, see Fig. 1 1(b) We schedule over the downlink and assume that 
each backhaul link is bidirectional and that both directions share the link capacity, but remark that all results can 
be readily extended to directional backhaul links. 


Definition 1 . User n is associated with up to two BSs: 

- The serving BS is denoted BS[n) and is defined as the BS that provides the highest SINR to user n. 

- The secondary BS is denoted BS(n) and is defined as the BS with highest SINR that has a backhaul link to 
BS(n) in C. 

Note that some users may not have a secondary BS. 

Packets destined for user n arrive at BS (n) and are stored in a queue Q n . Transmission for user n can be either 
single-transmitted by BS(n) or joint-transmitted by BS(n) and BS(n). For a packet to be joint-transmitted, it first 
has to be forwarded over the backhaul to BS(n), and stored in a queue Q n for joint transmissions maintained at 
both BS(ro) and BS(n). So a packet departs from Q n when it is single-transmitted or forwarded across the backhaul, 
and a packet departs from Q n when it is joint-transmitted. 


To illustrate this, consider the network in Fig. |2(a) with three users and three BSs. In Fig. 2(b) the primary and 
secondary BSs of each user are marked, with a solid and dashed line, respectively. Note that user 2’s secondary 
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BS2 user 2 



(a) example network 


Fig. 2. BS allocation and queues. 


~filQ2 



(b) primary (solid line) and secondary (c) the storage locations of the various 
(dashed line) BSs queues 


BS is BS1 and not BS3, although the latter is closer, since BS3 does not have a backhaul link to BS2. User 3 does 
not have a secondary BS because its primary BS does not have any backhaul connections. Fig. 2(c)| displays the 
various queues in play, and their locations. 

Both wireless packet transmissions and forwarding over the backhaul lasts exactly a single subframe. Throughout, 
we assume that a central processing unit determines the schedule for all BSs based on perfect knowledge on the 
network state. 

We consider a time-slotted model indexed by t, t = 0,1,..., where each time slot corresponds to a single 
subframe. Denote L n (t) and L n (t) the queue length of Q„ and Q n at the beginning of subframe t, respectively. 
Denote by W n {t) the number of new packets generated for user n at the beginning of subframe t. The W n {t) 
are assumed to be i.i.d. over time, independent between users, and have finite second moment. We denote by 


/in \t), and / Xn\t ) the number of packets transmitted towards user n in subframe t using single and joint 

transmission, and the number of packets forwarded across the backhaul, respectively. These are determined by the 
resource allocation at each subframe, see Section |III-A| for more details. The evolution of the queue lengths can 
then be written as 


L n {t+ 1) = L n (t ) + W n (t) - n\£\t ) - ^\t), (1) 

L n (t+ 1) = L n (t) + /4 3 >(f) - \t )• (2) 


A. Subframe model 

We now consider a single subframe consisting of scheduled blocks S = {l,...,jS} for each BS, see Fig. [3] 
and [v] we discuss how to allocate resources within a single subframe, which determines the 


In Sections 


IV 


j — 1, 2, 3. The evolution of the queue lengths (|T]» and ([2]> over time is then discussed in Sections [Vi] and VII 


A packet i which is single-transmitted requires scheduled blocks in a subframe of BS(n(z)), while a joint- 
transmitted packet requires scheduled blocks in the subframes of both BS(n(*)) and BS(n(i)). In the latter case. 
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Fig. 3. Example of a frame (corresponding to 1.4MHz LTE). 


the set of scheduled blocks used by BS(n(i)) and BS (n(i)) must have identical indices since JT requires both BSs 
to transmit on the same scheduled blocks. 

A packet i is characterized by the pair (n,/3), where n is the receiving user, (3 £ {0,1} indicates whether a 
packet is in Q n (f3 = 0) or in Q n ((3 = 1). The set of all packets is denoted I, \I\ = I. Given a packet i £ I, we 
denote by n(i), (3(i) its user type and queue type, respectively. 

When transmitted wirelessly, packet i is received successfully with probability p{i). Additionally, the success 
probability p(i) is independent of its allocated scheduled blocks, since we assume that the interference in all 
scheduled blocks is similar (due to frequency reuse 1). Note that if some scheduled blocks are unused, the 
interference is lower and better performance is obtained. We assume that p(i) is higher if (3(i) = 1 compared 
to (3(i) = 0, since when using joint transmission the two BSs configure their transmission parameters such that the 
signal combines constructively at the user, resulting in greater SINR when a packet is joint-transmitted. 

To simplify the presentation and due to space constraints, we make three assumptions: (i) forwarding a packet i 
over the backhaul is always successful; (ii) all packets are of same length and a packet transmission on a wireless 
channel requires one scheduled block from the subframe of each of its transmitting BSs; and (iii) the capacity 
of each backhaul link is K packets/subframe. In Appendix [B] we show that assumptions (i)-(iii) can be relaxed. 
Moreover, there we demonstrate the applicability of our algorithms to the case where a packet can be transmitted 
using one of several Modulation and Coding Schemes (MCSs). We also remark that all the results in this paper 
can be readily applied to the setting with infinite backhaul capacity, by setting K = S. In our simulation study 
(Section |VII|) we evaluate our algorithms for the case where multiple MCSs are supported. 


IV. The OFDMA Joint Scheduling (OJS) Problem 


We now formulate the joint scheduling problem. In order to describe the BSs involved in each transmission, we 
introduce 


{BS(n(*))| if (3(i) = 0, 

{BS(n(*)),BS(n(*))} if f3(i) = 1. 


If (3(i) = 0 then h returns only the serving BS, and if ti{i) = 1 it returns both the serving and secondary BS. 
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The function u : I x {0,1} h-» R + represents the utility of scheduling packet i over the backhaul (zt(i,0)) or 
wireless channel (u(i, 1)). Examples include throughput-based utility function ut and fairness-based utility function 
uf defined by 

■u T (i, 0 ) = 7 , u T (i, 1 ) = p(i), 

V*'/ 

u F (i, 0 ) = 7 , u F (i, 1 ) = logp(i), 


where 7 > 0 is a small constant that ensures packets are forwarded over the backhaul. Since here we consider a 
single-slot formulation, the utility of scheduling packets over the backhaul is not evident when the utility function 


is based only on wireless throughput; 7 compensates for this. In Sections VI and VII we use a queue-length based 
utility function 

UQ(i, 0) = max{L n( j) - L n(i) ,0}, 

u Q (i,l) = L n{i) p{i), 


(4) 


where L n and L n denote the queue length of Q n and Q n , respectively. Our model and analysis can in fact handle 
a wide range of utility functions such as those used in |27) . 

Based on the set of packets I and the utility function u, the centralized scheduler determines the set of wireless 
transmissions to take place in the upcoming subframe, and what packets to forward over the backhaul. 

The scheduler must also determine which scheduled blocks will be used for each packet transmission, such 
that for IT the scheduled blocks of the serving BS and secondary BS are aligned (i.e., have identical index). The 
scheduling decisions are represented using indicator variables z t £ {0, 1 }, yi £ {0, 1 }, and Xi S £ {0, 1 }, where 
z,j indicates if packet i is transmitted wirelessly, y, indicates if packet i is forwarded over the backhaul and x lH 
indicates if scheduled block s is used by packet i. The scheduler needs to solve the following integer programming 
problem (with 2 = y = (t/i) ie x, and x = (x is ) ie x jSe5 ). 

OFDMA Joint Scheduling (OJS) Problem: 

1 

max Ziu(i,l) + yiu(i,0) =: U{z,y ) 

x,y,z L ' 

i =1 


• Vi < 1? V i G X, 

(5) 

y; Zi<S,Wa£B ; ^ y t < K, Wl £ C, 

( 6 ) 

{i:a£h(i)} {i:h(i)=l} 

S 

y Xi S — Zi, V* £l; yi = 0, V« s.t. j3(i) = 1, 

S = 1 

(7) 

y x is <1, \/b £ B Vs £ S, 

( 8 ) 

{ i : b£h(i)} 

Zi £ {0,1}, Iti £ {0,1}, Vi £ X, 

(9) 

Xi S € {0,1}, Vi £ X Vs £ S. 

( 10 ) 


Constraint <[5]> ensures a packet is scheduled at most once and resides in a single queue; ([ 6 ]i ensures capacities in 
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(a) a wireless network with 3 users and 7 
packets (denoted 1..7) 


3 

2 


3 

6 


4 

5 
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(b) a possible allocation of 
the scheduled blocks 



(c) the resulting scheduled block graph 


Fig. 4. An example of a scheduled block allocation and the resulting scheduled block graph. 


each subframe and backhaul link are not exceeded; 0 ensures a scheduled block is allocated for each wireless 
transmission and packets in Q n cannot be forwarded over the backhaul; and ([8]) ensures that each scheduled block 
is used at most once in the subframe of each BS. 


To illustrate this problem, consider the network displayed in Fig. 4(a) This network comprises 3 BSs, 3 users and 
7 packets, numbered 1,... ,7. BS1 and BS2 are connected with a backhaul, so they can joint-transmit packets for 
user 1. Let S = 2, and assume that we want to allocate the scheduled block to achieve the following transmissions: 
(i) packet 3 should be joint-transmitted; (ii) packets 2,4,5 and 6 single-transmitted; and (iii) packet 1 forwarded 
over the backhaul from Q \ to Q\. This schedule can be obtained with the assignment displayed in Fig. 4(b) and 
the following solution to the OJS problem: 


2 = (0,1,1,1,1,1,0), y= (1,0,0,0,0,0,0), x = 

We now describe the complexity of the OJS problem. 


0 0 1 1 0 0 0 
0 10 0 110 


Proposition 1. OJS is strongly NP-hard even for instances in which all of the following hold: 

(a) u{i , 0) = 1, u(i, 1) = 1, V* S I; 

(b) L n < 1, L n < 1, Vn e Jf. 


The proof of Proposition [I] uses a reduction from the well-known problem of minimum edge coloring [ 191. The 
reduction demonstrates that even for cases with sufficient bandwidth to accommodate all packet transmissions in 
the BSs subframes, obtaining a feasible schedule where joint-transmissions use an identical set of scheduled blocks 
is equivalent to the well-known problem of minimum edge coloring G3- In Section [V] we use algorithms for 
minimum edge coloring when developing algorithms for OJS. 


V. OJS Problem - Algorithms 


In this section we develop algorithms to solve the OJS problem. First, we describe a framework for solving the 


OJS problem by decomposing it into problems related to knapsack and coloring, see Section V-A We then use 


this decomposition framework to develop efficient algorithms for OJS. In particular, in Section |V-B| we develop 



























TABLE I 

The performance and input required for different algorithms A Q js = [A JT k, A jtc ] 
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A JTK 

^4jtc 

Ratio 

Running time 

Input Gj 

JTK-MMK 

JTC-BIP 

a 

0(T M mk(I,B,C, S)) 

bipartite 

JTK-MAT 

JTC-BIP 

3A(G ,) 

0(CT mmk (I, 2,1, S)) 

any 

JTK-STA 

JTC-BIP 

A(Gj) 

0(T M mk(/.A(Gj) + 

1, A(G j), S)B 2 ) 

any 

JTK-PSP 

JTC-PSP 

a 

0{ 2 B + Tmmk(I , 

£ + 2 b ,C, S)) 

planar 
ser. paral. 


algorithms for instances where the joint transmission graph Gj (consisting of the BSs and the backhaul) is bipartite, 
in Section |V-C | we develop algorithms for instances where the joint transmission graph is planar and series-parallel. 


and in Section V-D we develop algorithms for general joint transmission graphs. Note that joint transmission graphs 
encountered in practice need not always be bipartite or planar and series-parallel. However, if this is the case, by 
using an algorithm that exploits the structure of the graph, we can guarantee lower computational complexity and 
more accurate results. The algorithms for the general case in Section |V-D| are based on those for bipartite graphs 
in Section IV-BI 

We denote the approximation ratio of a given algorithm by a (0 < a < 1). If the algorithm is optimal, we have 
a = 1 . 


A. Decomposition Framework 

From Proposition |T] we see that, unless P=NP, an efficient optimal algorithm for general instances of OJS does 
not exist. In order to develop efficient approximations for the general case and optimal solutions for a subset of 
instances, we present two additional scheduling problems and explore their relation to OJS. These two problems 
are obtained by partitioning OJS into two parts, exploiting the fact that x only appears in 0, ([8]), and ( [T()| . 

Joint Transmission Knapsack (JTK) Problem: 

max U(z,y) 
z.y 

s.t. 0,0,0, 3x: 0. 0, ([10} hold 
Joint Transmission Coloring (JTC) Problem: 

given z.y , find x s.t. 0. 0, ( fTO} hold 

Note that the JTK problem resembles an assignment problem rather than a knapsack problem. We remark that 
this is due to‘ the assumptions that all packets have the same length and use a single MCS. We show in Appendix [B] 
that when relaxing these assumptions, JTK indeed generalizes to a knapsack-like problem. 

We use A jtk and /1 JTC to denote an algorithm for JTK and JTC, respectively. A specific algorithm for problem 
P is denoted P-D where D identifies the algorithm. For instance, we write Ajtk = JTK-GREEDY if we solve the 
JTK problem using a greedy algorithm. An instance of problem P consists of specific values for all variables in its 
constraints, except for x, y, and z (JTK) and x (JTC). 
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Fig. 5. The decomposition framework leading to algorithm Tqjs =1 /I.itk- Aitc I- 


The JTK problem differs from OJS in that it does not attempt to find x but guarantees that such x exists for 
its solution z, y. The JTC problem then finds x given z and y, which is later shown to correspond to a coloring 
problem. It is ensured that if a solution exists for JTK, the corresponding JTC problem instance can also be solved. 
Thus we can decompose OJS by first solving JTK to obtain z,y and then solving JTC to obtain x, see Fig. [5] 
(some details in the figure are explained later). In Sections V-B||V-D we identify instances where the existence of 
a feasible x is guaranteed without the need to find x, and use this to develop efficient algorithms for OJS. 

If algorithms Ajtk and /4 n c are used in this decomposition to solve JTK and JTC, respectively, we denote 
the corresponding OJS algorithms as ,4 0 | S =[.4 | T k, A nc ]. As we shall demonstrate in this section, making this 
decomposition allows us to find efficient algorithms for solving OJS. Throughout the paper only optimal Ajtc 
algorithms are considered. The following lemma immediately follows. 

Lemma 1. If Ajjk is an a-approximation algorithm for the JTK problem and Ajtc an optimal algorithm for the 
JTC problem, the algorithm Aqjs —[Ajtk, Ajtc] is an a-approximation for the OJS problem. 


We introduce the following definitions that will be used to solve JTC. Recall from Section III that Gj = (B. C) 
denotes the joint transmission graph that describes what BSs are connected by a backhaul link. 

Definition 2. The scheduled blocks graph Gsb = Gsb(G j, z) = ( Vsb,Esb) of a JTC instance is defined by 

V SB = Vj U {B + 1,..., 2 B} 

and 

E$b = |{BS'(n(f)),BS(n(i)) + B} | z* = 1 and /3(i) = oj U j {BS(n(i)), BS(n(i))} | z* = 1 and j3{i) = lj. 

To interpret the definition of Gsb- recall that z* denotes whether packet i does a wireless transmission (z, = 1) 
or not (zj = 0). So each edge in the scheduled block graph represents a wireless transmission, connecting the BSs 
that are engaged in the transmission. To avoid self-loops, we introduce dummy vertices B + 1,..., 2 B in case these 
BSs are involved in a single-transmission. Note that while Gj is a simple graph, Gsb need not be. As an example 


for Gj and Gsb consider the network depicted in Fig. 4(a) Using the solution described in the example at the 


end of Section [TV] the resulting scheduled block graph is depicted in Fig. 4(c) 


We now show that, using the scheduled block graph, JTC can be rewritten as an instance of the well-known 
edge-coloring (EC) problem |10|, and solved accordingly. The input to the EC problem is a graph G = (V, E) and 
the output is a coloring on the edges that uses a minimum number of colors. 
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Lemma 2. JTC is equivalent to finding an edge coloring using at most S colors on Gsb = (Vsb,E sb ). 

As a consequence, JTC can be solved optimally by invoking an optimal algorithm /l H c for the EC problem on 
Gsb- Note that the EC problem is NP-hard. 

Table [I] summarizes the different options we will describe in sections |V-B||V-Dl for solving OJS using the 
decomposition framework. 


B. Algorithms for Bipartite Network Graphs 

We now develop algorithms for OJS instances in which Gj is bipartite. The results in this section can be used to 
solve such networks, and will provide the building blocks for the algorithms for general joint transmission graphs 


in Section V-D We start by describing an algorithm for JTK and an algorithm for JTC, and show how using them 
in the decomposition framework will result in an approximation algorithm for OJS. We require the following two 
lemmas. Denote by A(G) the maximum vertex degree of G 


Lemma 3. If Gj is bipartite, then Gsb(Gj,z) is bipartite for every z. 

Lemma 4. If Gsb A bipartite and A {Gsb) < S, then 3x such that (f7]> , ([8|>, (flQ| hold. 

We now describe an algorithm for JTK based on the well-known Multidimensional Multiple-choice Knapsack 
(MMK) Problem pT) . Observing the formulation of MMK in |2 1 1 , the input to MMK is a subset of the input to 
OJS. We define the algorithm Ajtk =JTK-MMK as simply running some algorithm AyiviK for solving MMK with 
a running time of Tmmk (for different Ammk algorithms see Table 0 and show that it solves JTK for bipartite 
networks. 


Lemma 5. If Gj is bipartite and algorithm A MMK is an a-approximation algorithm for MMK, A JTK = A MMK is 
an a-approximation algorithm for JTK. 

Next, we describe an algorithm for JTC when the network graph is bipartite, by exploiting the connection to graph¬ 
coloring problems from Lemma [2] Let JTC-BIP be the edge coloring algorithm from [TOj. Using Lemma [3] Gsb is 
bipartite and since also A(Gsb) < S it follows from 1101 that JTC-BIP finds an edge coloring using at most S colors. 
Using Lemma [2] we conclude that JTC-BIP solves JTC. The running time of JTC-BIP is Od-EsBl log A(Gsb)) = 
0(BS log S). The following theorem is the main result of this section. 

Theorem 1. For bipartite networks, if Ajtk =JTK-MMK is an a-approximation for JTK, then Aojs =[JTK-MMK, 
JTC-BIP] is an a-approximation for OJS. 


C. Algorithm for Planar Series-Parallel Graphs 

We now develop an algorithm for OJS instances in which Gj is planar and series-parallel. We describe algorithms 
for JTK and JTC, and use them in the decomposition framework to devise an approximation algorithm for OJS. In 
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TABLE II 

Algorithms for MMK 


^ 4 .mmk 

Ratio 

T M mk(I,B,C,S) 

DP (21 

Optimal 

0(S ib+c;) /(B + C)) 

PTAS 1261 

1/(1 + 4 

o(/(( B+c ^/ ) 

Greedy I 2 I | 

00 

O(llog(l)) 


this section we use similar ideas to those in Section V-B We need the following lemma whose proof is similar to 
that of Lemma [2 


Lemma 6. If Gj is planar and series-parallel, Gsb is planar and series-parallel for every z. 

We first describe Algorithm JTK-PSP that solves JTK when Gj is planar and series-parallel. The algorithm 
uses .4 mmk to solve an MMK instance defined as follows. The number of dimensions is D' = D + |i5 0 dd| where 
Bodd = {& C B : \B'\ is odd and > 3}. The capacity for each new dimension associated with a set B' £ B 0 && is 
iS(|£? , | — l)/2. The weight in each new dimension for each (i,r) such that r > 0 and h(i,r) C B’ is set to T(i,r); 
for all other cases it is set to zero. The algorithm concludes by scheduling packets for transmission according to 
the configurations selected in the solution returned by ,4 mm k- 

We note that |B 0 dd| = 0(2 B ) and therefore in general Algorithm JTK-PSP may be impractical due to a very 
large running time. This algorithm is therefore more appropriate for small B. We now show that in some instances 
the running time can be improved. Note that if a set B' £ B,m has no more than \B'\ — 1 edges in C that connects 
two nodes in B' , it can be removed from £> 0 dd; if the number of such sets is large this can significantly decrease 
the running time. 

For networks that are planar and series-parallel, JTC can be solved using the edge-coloring algorithm from f36) . 
We call this algorithm JTC-PSP, and note that its running time is 0( BA(Gsn)) = 0(B ■ S). The following theorem 
applies the decomposition framework to planar and series-parallel networks. 

Theorem 2. For planar and series-parallel networks, if JTK-MMK is an a-approximation for MMK then [JTK-PSP, 
JTC-PSP] is an a-approximation for OJS. 


D. Algorithms for General Graphs 

We now develop algorithms for general OJS instances, without imposing any conditions on Gj. We start with 
describing two approximation algorithms for JTK. For each approximation algorithm we show how using it in the 
decomposition framework will result in an approximation algorithm for OJS. 

First, we describe Algorithm JTK-MAT which is based on computing a matching. For each {a, b} £ C, the 
algorithm solves an instance of JTK defined by a network that has only two BSs a,b and the backhaul link with 
capacity K. Only packets that can be scheduled in such network are considered, and .4 mmk with B = 2 and (7=1 
is used to solve this limited instance. Each edge in C is assigned a weight equal to the total utility obtained when 
solving its limited JTK instance. Then, maximum weighted matching is found and the union of all solutions for 
edges in the matching is returned. This solution is feasible for the general JTK problem. 
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Algorithm JTK-MAT Based on matching 
l: for b € B do 

2: If b has no backhaul link, run Ammk to solve a JTK instance with B' = {6}, I' = {i £ I : h(i) = B'}, 

C' = {} and set z, and yi as determined by A M mk- 
3: for e = {a, b} G C do 

4: Run A M mk to solve a JTK instance with B' = {a, 6}, I' = {t £ 1: h(i) C B'}, C = {e} 

5: Assign U(z,y) found by A M mk as a weight for e 

6: Compute maximum weight matching on Gj and store the result in the edge set 8 

1: Set Zi = 1 (y, = 1) only if 3e G 8 such that z, = 1 (yi = 1) in the solution returned in line [4] for edge e 


Theorem 3. If algorithm Ammk I s an ca-approximation algorithm, then Aqjs =[JTK-MAT (Ammk), JTC-BIP] 
is a (2a) / (.3 A (G j)) -approximation algorithm for OJS. 

The maximum weight matching algorithm from |l6| that takes 0(|£7||T^| log |Vj) time can be used in line [6} in 
this case the running time of Algorithm JTK-MAT is dominated by Ammk in line [4] Therefore, the running time 
of Algorithm JTK-MAT is 0(CTmmk(I, 2, 1, S)), where Tmmk is the running time of Ammk (see Table |5). 

We now describe Algorithm JTK-STA which iterates over star subgraphs. It is similar to Algorithm JTK-MAT, 
but iterates over the vertices b £ B instead of the edges e G C. The approximation ratio of JTK-STA is better than 
that of JTK-MAT, but its running time is worse. For each vertex b, a JTK instance is constructed using only b and 
its neighbors in Gj. The solution of this instance, z,y, is assigned to b. Next, the algorithm finds the vertex 6 max 
associated with maximum total utility. If 2max■ 2/max is the solution associated with 6 max , the algorithm schedules 
the packets indicated by z max , 2/max- The vertex 6 max and its neighbors are removed from B. This process is repeated 
until B is empty. 

Note that after the first vertex is removed from B, in order to update the weights it is sufficient to consider 2-hop 
neighbors of & max in line [14] since weights of other vertices remain unchanged. The running time of JTK-STA is 

0(Tmmk(J, A (Gj) + 1, A (Gj),S)B 2 ). 

The following theorem proves that OJS can be solved approximately when JTK-STA is used in the decomposition 
framework. 

Theorem 4. If algorithm A MMK is an a-approximation algorithm, then, Aqjs =[JTK-STA (Ammk), JTC-BIP] is 
an (a/ A(G j))-approximation algorithm for OJS. 


VI. Queueing Dynamics 

The OJS problem discussed in Sections IV and [V] schedules transmissions within a single subframe. We now 
expand the scope to multiple subframes, where packets arrive and depart over time, and study the evolution of the 
users’ queues. Our objective is to identify a scheduling policy that is maximum stable (or throughput optimal). 
Under such a policy, the queue-length process is positive recurrent for any arrival for which a stabilizing policy 
exists, see, e.g., (13), (25), (30). We prove that by using a specific utility function and an algorithm for solving 
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Algorithm JTK-STA Based on star subgraphs 
l: function SOL-STAR(6, B', C, X') “ ~ 

2: Run Ammk to solve a JTK instance defined with B' = {6} U {a : {a,b} £ C'}, C = {{a, 6} : {a, 6} £ 

C', a £ B', b£ B'}, T = {i£l' : h(i) C &'} 
return z,y as determined by Ammk 

3: for b £ B do 

4: Z[b],Y[b\ <- SOL-STAR(L B, C, X) 

5: Assign U(Z[b],Y[b]) as a weight for b in B 

6: Initialize B" <- B; C" 4- C; X" £-X 

i: repeat 

B: Find the vertex 6 max in B" with maximum weight 

9: X 4- {i £ I" : 2[/) m ax]i + ^ [&max]i = 1} 

10: for all i £ X, 4 Z[b mSLX ^i, y r 4 y [6 ln;i x]? 

11: J £- {a e B" : 36 £ B", {& max , b} £ C", { a , b} £ C"} 

12 : X" ^1"\1-B" ^B"\{a: {a,6 max } G C") 

13: Remove from C" edges with an endpoint not in B" 

14: for b £ J do 

15: Z[b\,Y[b\ 4- SOL-STAR16, B", C ", X") 

16: Update U(Z[b\,Y[b}) as a weight for b in B" 

17: until B" is empty 

18: return z" 


the OJS problem (we refer to this combination as a scheduling policy), we obtain a MaxWeight-like scheduling 
policy (see, e.g., ®), which is throughput optimal. 

Let L{t) = (Li(f), Li(t ),..., L n (t), L n (t)) denote the queue length process at time t, so {L(t )} t >o is the 
stochastic process that tracks the queue-length evolution over time. We solve the OJS problem in each subframe, 
given a certain utility function. Let us denote by y(t) = (t/i(f))igi and z(t) = (zi(t))i^x the solution of OJS in 
slot t. Here ( y(t),z(t )) can represent both an exact solution or an approximation. 

Recall that W n (t) is the number of packets generated for user n at time t. Let A„ = E{W„(0)}, and define 
A = (Ai,...,Ajv) to be the arrival rates. We denote by fin\t;z(t)), yn\t;z(t)), and ^n\t\y(t)) the number of 
packets transmitted towards user n in subframe t using single and joint transmission, and the number of packets 
forwarded across the backhaul, respectively, given solution ( y(t),z(t )) of OJS. 

Denote X n {t) the set of packets in Q n at time t (so \X n (t)\ = L n (t )), and X n (t) the set of packets in Q n . The 
fin , j = 1,2,3 can be written as 

Vn ) ( t '> z )= (f; y) = 

where the Y it ~ Ber(/v(<)) are mutually independent Bernoulli distributed random variables that represent whether 
packet transmissions are successful. For notational convenience, we write yn\t) to represent the transmission rates 
at time t. 

It is readily seen that the joint queue-length process {X(f)} t >o is Markovian. 

We now analyze the traffic intensity that can be sustained by the queueing system described by o and (|2j. The 
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Fig. 6. Simulated network topologies. |(a)| SINR for the cluster of 3 BSs. The dark 

blue area denotes the location of inter-cell users. |(b)|(c)| Two 7 BS topologies. The big. 7. SINR for the cycle topology for (a) single¬ 
color of the nodes illustrate the throughput gains of users for a sample run. Red transmissions and (b) joint-transmissions, 

users benefit from JT and the throughput of yellow and black users remains roughly 
the same. Black users have throughput of 0 both with and without JT. 


stability region of a particular policy is defined as the set of all arrival rates such that the {L(f)} t > 0 is positive 
recurrent is called the stability region of this particular policy. The capacity region of a network is defined as the 
union of the stability region over all policies. If the stability region of an algorithm OJS-ALG and utility function 
u is equal to the capacity region, we say that policy (OJS-ALG,u) is throughput-optimal. 

In order to investigate the network capacity region in more detail, we first introduce some definitions. We denote 
by fn 1 ' ) the rate (long-term average traffic flow) of single-transmission packets for user n, by /^ the rate of joint- 
transmission traffic for user n, and by /,, the rate of user-n traffic sent across the backhaul. Define the vector 
/ = f[ 2 \ f[ 3 \ • • •, fx\ /,y ■*■ /jv^)- Then, for a given arrival rate vector A, the set of all A-admissible traffic 

flows can be defined as 


F x = 


{/ e Ry V 



(ll) 


That is, a flow is A-admissible if for all queues Q n , Q n , n = 1..., N, the traffic arrival rate is equal to the departure 
rate. 

We now introduce the set of all arrival rate vectors such that at least one A-admissible flow can be realized: 


A = {a G R 3 + n I 3f G F x 3r G conv(f?,)/« < r« if /« > 0, n G AT, j = 1,2, 3 }, 

where conv(i?) denotes the convex hull of R, the set of all rates across the various links that can be achieved in 
saturation. 

We now show that any A G A can be stabilized, and that any A outside of the closure of A cannot. For stabilizing 
A G A we use the policy (OJS-OPTuq), where OJS-OPT represents any algorithm that solves OJS exactly, and 
uq is the queue-length based utility function from 0 - The following theorem then implies that A is indeed the 
capacity region, and that (OJS-OPT,«q) is throughput-optimal. The proof relies on a standard drift argument using 
a quadratic Lyapunov function (30J. 

Theorem 5. Let A G A, then the queue-length process is stable under policy (OJS-OPT, uq). If A A, then there 
exists no policy that stabilizes the network. 


SINR (dB) 
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The reason for choosing the queue-length based utility function |4} becomes clear when substituting it into the 
objective function of OJS. 


N 

U(z) = ]T (L„(i)E{^)(f)} + L n (t)E{^\t)} + CL n (t ) - L n {t))p^{t)). 


( 12 ) 


n—1 


This yields the objective function of the celebrated MaxWeight scheduling algorithm. This algorithm, first introduced 
in pO] , has been shown to be throughput-optimal in a wide range of settings, see, e.g., Note that although 

our objective function of maximizing the queue-weighted throughput is similar to that used in traditional work on 
MaxWeight scheduling, the constraints of the OJS problem are markedly different. Specifically, the MaxWeight 
scheduling literature is typically concerned with maximum (weighted) set problems, which are fundamentally 
different from the OJS problem. 

Since the maximization for utility function 0 is a specific instance of the OJS problem, it follows from 
Proposition [I] that solving this problem is NP-hard. Thus, using an optimal algorithm OJS-OPT (as in Theorem [5]) 


is typically not feasible in practice for general graphs. In Section VII we investigate the performance of a wider 
set of (suboptimal) algorithms via simulation. Theoretical results on the capacity loss for general algorithms will 
be the subject of future work. 


VII. Numerical Results 

We conducted a simulation study to evaluate the performance of the various algorithms introduced in Section |V] 
Throughout this section we consider the case where a packet can be transmitted using one of several Modulation and 
Coding Schemes (MCSs); Details on extending the algorithms to support several MCSs are given in Appendix [b] 
The simulation results provide insights on the network-level benefits and tradeoffs of JT under various network 
scenarios. 


A. Simulation Setup 

OJS algorithm. We implemented the four algorithms presented in Table [T] The majority of the algorithms 
and the queueing dynamics are implemented in Python, while the ,4 MMK procedure is written in C. We did not 
implement any JTC algorithms, since it was proved in Section [V] that JTK guarantees that a feasible solution always 
exists. Furthermore, as we show in Section VII-B1 the greedy algorithm shown in Table [II] performs well in most 
considered scenarios. 


Network setup. We consider three network topologies. In Section VII-B2 we analyze a network of 3 BSs, with 
backhaul links between each pair (Fig. |6(a)| i. In Section |VII-B3[ we look at two different backhaul topologies for 
a 7 BS network, shown in Fig. 6(b)||6(c) The backhaul capacity (BC) is the same for all links. 

We use a fixed packet size of 73 bytes, and the backhaul capacities are normalized to units of packets/subframe. 


The distance between neighboring BSs is 700 m. The height of each BS’s antenna is 20 m. The BSs’ transmit 
power is 39 dBm and 30 dBm for the 3 and 7 BSs network, respectively. Lower transmission power is used for 
the larger network, since more BSs transmit interfering signals. 
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We simulate N = 20 users for the 3-cluster, and N = 50 for the 7-cluster topologies. The users are placed 
uniformly within a circle that contains the entire simulation area, and with radius 1050 m. 

Wireless model. We set S = 50 scheduled blocks, corresponding to a 10MHz LTE system. Once the location for 


a user is determined, the received power level from each BS is computed based on the Hata propagation model 1181 
which was shown suitable for LTE in urban areas (29) . The power levels from the different BSs are used to compute 
the SINR for single and joint transmission to the user. The SINR values for single and joint transmissions for the 
cycle topology are plotted in Fig. [7] Given the SINR values, the success probability for single and joint transmission 
is computed for each MCS (QPSK-1/2, QAM64-1/2, and QAM64-3/4) using data taken from (4). 


Queueing Dynamics. The queueing dynamics are implemented as in Section VI Unless otherwise noted, packet 
arrivals for users follow a binomial distribution with n = 3 and p = 0.5. In every subframe, an algorithm for solving 
OJS is executed with the utility function <|4}. Throughout the simulation, we track the normalized throughput of 
a user, defined as the fraction of arrived packets that have been successfully transmitted. The average normalized 
throughput is computed over all users. The simulation duration is 1,000 subframes and each data point is obtained 
by averaging over 1,000 runs. 

We also distinguish the performance of inter -cell users from that of the infra-cell users in order to evaluate the 
benefits of CoMP JT for both types. In this section, we refer to inter-cell users as those whose power levels from 
two BSs is above a threshold. This threshold is determined numerically by observing the physical location of these 
users (illustrated by the darker regions of Fig. |6(a)|>. 


B. Simulation Results 

1) Performance of the Approximation Algorithms: In practice, approximation algorithms may perform signifi¬ 
cantly better than their guaranteed approximation ratios. Hence, we carry out a single-subframe evaluation with 
the goal of isolating the performance of the algorithms from the long-term effects of the queueing dynamics. We 
consider the OJS formulation with the throughput utility function (3). We consider all algorithms, for two different 
topologies with 3 BSs: complete graph and a bipartite graph. 

When an optimal Ammk subroutine is used, JTK-MMK and JTK-PSP are optimal under the complete graph and 
bipartite topologies, respectively, by Theorems |T| and [2] These two algorithms, are therefore, used as benchmarks. 
JTK-MAT and JTK-STA operate on general networks, but only consider a subset of the backhaul links. Therefore, 
they achieve a fraction of the optimal utility, denoted utility ratio in Fig. [8] We vary the number of users from 1 
to 80, reflecting the range of users that can be expected in small cell deployments. In each run, a set of items X is 
sampled randomly. To obtain a single point, 10,000 iterations are averaged. 

We first use the optimal DP algorithm for A MMK , for which a = 1 (Table [II]). Since the maximum vertex degree 
is A (Gj) = 2 for both topologies under consideration, JTK-MAT and JTK-STA are 1/3- and 1/2-approximations, 
respectively (Table [I). For the complete graph topology, JTK-MAT achieves a utility ratio of 0.6 at its worst, while 
JTK-STA does better with a ratio of 0.8 (Fig. |8ja)). Similar insights hold for the bipartite topology in Fig. [8j b). 
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— JTK-STA 

+--+ G-JTK-STA 

* . G-JTK-PSP 

—• JTK-MAT 
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Number of users Number of users 

(a) complete graph (b) bipartite graph 


Fig. 8. The ratio between the optimal utility for the OJS problem for a single subframe, and the utility obtained by the different algorithms 
under the two topologies. 



Backhaul capacity 
(a) inter-cell 


Backhaul capacity 
(b) intra-cell 


Fig. 9. The throughput of both intercell and intracell users for a range of backhaul capacity levels. 


In addition to using the optimal Ammk, we ran the same simulations for the case when a greedy algorithm is 
used for Ammk (Table |HJ, denoted with the prefix “G-” in Fig. [8] In this case the approximation ratios no longer 
hold, as there are no performance guarantees for the greedy algorithm. However, we found that when 4 mmk is 
solved greedily, the algorithms are very close to optimal for more than 10 users. Moreover, the running time of the 
greedy algorithm is significantly lower than the duration of a subframe. Due to their improved running time, we 
only use the greedy version in the following sections. For clarity, we omit the “G-” prefix. 

2) Impact of Backhaul Capacity: Backhaul links are typically expensive to deploy, and operators frequently have 
to lease them. Therefore, it is important for the operator to strike a balance between improving performance and 
containing backhaul costs. To obtain a better understanding of the required backhaul capacity, we evaluated its 
impact on the long-term throughput of the the queueing system. 

In Fig. [9] the user-averaged normalized throughput is shown when the backhaul capacity between each pair of 
BSs is scaled from 0 to 6 packets/subframe. Inter-cell users in particular gain from JT (Fig. [9] a)) when network-level 
behavior is considered. A 28% throughput gain is observed for those users with the addition of backhaul. Half of 
this gain is achievable with 1 unit of backhaul capacity, while 2 units realizes 80% of the potential gains. However, 
intra-cell users gain 5% throughput. As cell sizes become smaller, the portion of inter-cell users increases and the 
overall gain from using CoMP JT will be higher. 

The achieved throughput depends on the used algorithms. The largest benefits are possible with the JTK-PSP 
and JTK-STA algorithms, which utilize 3 or 2 backhaul links in every subframe, respectively. It is also observed 
that in clusters of this size, JTK-STA performs as well as the optimal JTK-PSP, despite its lower running time. 
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JTK-MAT requires more backhaul capacity to achieve the same throughput. Overall, in this case JTK-STA is the 
best choice while JTK-MAT can be used to reduce computational resources at the cost of a larger investment in 
infrastructure. 


3) Impact of Topologies: We now consider the star and cycle topologies with 7 BSs, illustrated in Fig. 6(b) 


and 6(c) respectively. Since both of these topologies are bipartite, we use the optimal algorithm JTK-MMK in this 
section. 

To study the impact of network topology on JT, we study the throughput gains that are obtained when backhaul 
links with a 3 unit capacity are introduced in the different topologies. We invoke the algorithm for the same user 
placement and average the results. The throughput gains of individual users under the cycle and star topologies are 


illustrated in Fig. 6(b) and 6(c) respectively. Overall, we see that in each topology there are 4 users that observed 
an increased throughput from JT while the throughput of the other users remained very similar. While we do not 
list here the change in throughput for each individual user, the trend is as follows. The red users have throughput 
of 0 without JT, since the SINR for a single transmission to these users leads to a packet transmission success 
probability of zero. With JT, the throughput of the red users is very close to 1 and the throughput of each yellow 
users is reduced by at most 0.001. 

In both topologies we see improvement in the throughput of 4 users alongside a negligible decrease in throughput 
of some of other users. However, we observe that the users that benefit from JT (i.e., red users) are different for the 
different topologies. The throughput of users B,C,D, and E increased in the cycle topology but for the star topology 
the throughput of users B,C,F, and G increased. Clearly, this is to be expected as the users who gain from JT are 
those that reside between two BSs that are connected via a backhaul link. Additionally, we notice that in each 
topology there are two users that exhibit a throughput of 0 both with and without using JT. For the cycle topology 
(Fig. |6(b)| the users are F and G, while for the star topology (Fig. 6(c)| ) the users are D and E. To conclude, the two 
different topologies result with different distribution of the throughput among the users that are located between 
the BSs. Since the users are located somewhat evenly in the simulation area, the number of red users is similar for 
both topologies. However, for a different user placement a specific topology may benefit a larger number of users. 

The impact of the topology on the stability region of user A from Fig. 6(b)|6(c) is illustrated in Fig. 10(a) 


The aggregate queue size at the end of the simulation run, under different arrival rates and backhaul capacities, 
highlights that this user’s queues can be stabilized for higher arrival rates under the star topology. This behavior is 
representative of other inter-cell users. 

To further study the performance of our algorithms under the different topologies we plot in Fig. 1 10(b) and 


Fig. 10(c) the throughput as a function of the arrival rate, for the cycle and star topologies, respectively. In both 
figures JT improves the network throughput by 9%, even for low loads. This can be explained by the observation that 


certain inter-cell users may never receive a packet through single-transmission. For the cycle topology (Fig. 10(b) i, 
the performance of JTK-STA and JTK-MAT is comparable, despite the different approximation ratio. This is 
because in the cycle topology a matching may include 5 out of the 6 available backhaul links, making JTK-MAT 
comparable to JTK-STA. It turns out that for higher rate values, JTK-MAT performs slightly better than JTK-STA. 
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Arrival Rate (usr pkts/subframe) Arrival Rate (usr pkts/subframe) 



Arrival Rate (usr pkts/subframe) 


(a) a user’s queue sizes 


(b) Throughput for cycle toplogy (c) Throughput for star topology 


Fig. 10. |(a)| The final sum of queue sizes of user A from Fig. |6(b)| as the arrival rate is varied. Larger backhaul capacities (BCs) keep the 
queues bounded for higher arrival rates; [(b)|(c)| Throughput obtained by the proposed scheduling algorithms vs. the arrival rate for the star and 
cycle topologies. 


Therefore, an operator may choose to run both JTK-MAT and JTK-STA and select the schedule that yields the 
highest throughput in each subframe. Studying such algorithm is out of scope and left for future research. In the 


star topology (Fig. 10(c) i, since JTK-STA is an optimal algorithm, we omit the curve for JTK-MMK from the 
figure. For lower arrival rate the performance of JTK-MAT is very close to optimal; As the arrival rate increases 
JTK-STA is clearly favorable. This is expected since for the star topology the matching computed by JTK-MAT 
will include at most one backhaul link. 

The joint impact of backhaul capacity for star and cycle topologies is shown in Fig. 1 11 (a) | and Fig. 1 11 (b)| 


respectively. Recall that JTK-STA is optimal in Fig. 11(b) and note that, compared to the optimal JTK-MMK in 
Fig. 1 l(a)| bakchaul capacity of 3 or more is beneficial only for the cycle topology. This is due to the fact that 
users residing between the center BS and another BS exhibit greater interference than users that reside between 
two non-center BSs. Therefore, in the star topology, such users gain significant increase in their throughput by just 


forwarding their packets. As we observed in Fig. 6(c) there are two such users and since the arrival rate is 1.5, a 
backhaul capacity of 3 suffices for the star topology. However, in the cycle topology, while mainly two users gain 
a significant throughput increase, other users which can get their packets using single transmission can still benefit 
from JT. Therefore, increasing the backhaul capacity further improves the overall throughput. 

Fig. 11(a) and Fig. 11 (b)| also demonstrate the relationship between the network topology and the performance of 
our algorithms. For the cycle topology (Fig. 1 11 (a)[ >, a matching may include 5 out of the 6 backhaul links, for this 
reason JTK-MAT is relatively close to optimal and even surpasses JTK-STA for most backhaul capacity values. On 
the other hand, for star topology (Fig. 11 (b)| >, while JTK-STA is guaranteed to be optimal, the matching computed 
in JTK-MAT contains at most one backhaul link and therefore the algorithm obtains low performance. Note that 
JTK-MAT with backhaul capacity 6 performs closely to JTK-STA with backhaul capacity 1 since it utilizes at most 
one out of the 6 backhaul links in every subframe. 


VIII. Conclusions 

In this paper, we considered a cellular network with Coordinated Multi-Point (CoMP) Joint Transmission (JT) 
capabilities that allow multiple BSs to transmit simultaneously to each user. We first formulated the OFDMA 
Joint Scheduling (OJS) problem of determining a subframe schedule and deciding if to use JT. By exploiting the 
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(a) Throughput for cycle toplogy 

Fig. 11. Throughput obtained by the proposed scheduling algorithms 
rate of 1.5 packets per user per subframe. 



(b) Throughput for star topology 

vs. the backhaul capacity for the star and cycle topologies, for an arrival 


characteristics of this problem, we developed efficient scheduling algorithms for bipartite graphs and approximation 
algorithms for general graphs. 

We then considered a queueing model that evolves over time. In this model, we proved that solving the OJS 
problem with a specific queue-based utility function (in every subframe) achieves maximum throughput in CoMP- 
enabled networks. 

Via extensive simulations we showed that the bulk of the gains from CoMP with JT can be achieved with low 
capacity backhaul links. We showed that our algorithms distribute the network resources evenly, increasing the 
inter-cell users’ throughput at only a slight cost to the intra-cell users. 

This paper is the first step towards a rigorous, network-level understanding of the impact of cross-layer scheduling 
algorithms on CoMP networks with JT. In future research, we will extend the model by allowing more than two 
BSs to jointly transmit, and by allowing longer backhaul delays. Moreover, we will apply our techniques to CoMP 
related technologies such as network-MIMO, multi-cell MIMO, and MU-MIMO. Finally, we will study the design 
considerations of the backhaul network and the impact of decentralization on the performance. 
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Appendix A 
Proofs 


(Proposition [T]i To prove Proposition [I] we use the following definition: 


Definition 3. The chromatic index of a graph G [19], x'(G), is the number of colors required to color the edges 
of G such that no two adjacent edges have the same color. 


It is known by Vizing’s theorem that for every simple graph G, yj{G) = A(G) or }f(G) = A(G) + 1, where 
A(G) is the maximum vertex degree of G. The Minimum Edge Coloring Problem (MECP) 1191 is to determine 
whether x'(G) = A(G) or x'(G) = A(G) +1. It is well-known that MECP is NP-hard |19|, therefore to complete 
the proof we present a polynomial-time reduction from MECP to OJS. 

Given a simple graph G = ( V , E) with maximum vertex degree A(G), we now describe how to construct an OJS 
instance. We set B = V. For each edge {ui, V 2 } € E, we add a user n to J\f for which BS(n)=vi, BS(ro) = V 2 , and 
there exists only a single pending packet in Q n . Thus, for every packet i we have /3(f) = 1. The utility is defined 
as u(i, 1) = 1 and u(i, 0) = 1 for every i. Also, we set S = A(G), K = S , and C = {{a, b}, a,b £ B A a f b}. 
Note that for the constructed instance (a)-(b) stated in the lemma hold. 

We now show that the optimal solution to the OJS instance has a utility of \E\ if and only if yf(G) = A(G). 
Let x*, y* , z* be an optimal solution to the OJS instance with a total utility of \E\. Due to the utility u used and 
since the total number of packets is \E\, z* = 1, Vi Consider an edge e = {^ 1 ,^ 2 } G E and a pending packet i' 
such that h(i') = {v\,V 2 \- Due to there is exactly one s' for which Xi> s i = 1. We assign e the color s' and 
continue the process for the remaining edges and packets. Since 1 < s' < S = A(G), at most A(G) colors are 
used. Since ([8]) also holds, no two adjacent edges are colored using the same color s. We showed an edge coloring 
with at most A(G) colors, xj(G) = A(G). The other direction, namely showing that if x'(G) = A(G) then the 
optimal solution to the OJS instance has a utility of | E\, can be proved similarly. ■ 

(Lemma [lj: The proof immediately follows from the definitions of the JTK and JTC problems. ■ 

(Lemma [ 2 }: Given a solution x' to a JTC instance, we now define a coloring on Gsb that uses at most S 
colors. Observe that by Definition [2] and since constraint 0 holds, there exists a one-to-one mapping from every 
pair (i,s) such that x' ts = 1 into an edge in Gsb- This mapping defines an edge coloring using at most S colors 
(1 < s < S). Since constraint 0 holds, no two edges of the same color touch a vertex in Gsb- The other direction, 
namely, finding a solution x' to JTC, given an edge coloring on Gsb- can be proved similarly. ■ 

(Lemma [3|: Using Definition Q, it is clear that the subgraph G' = (V',E') of Gsb defined by V = B and 
E' = {{a, b} £ Esb '■ a, b £ B} is bipartite. Since each b £ Vsb \ V' has at most one neighbor, Gsb is bipartite. ■ 
(Lemma[4j: Using the result from jlO] and since Gsb is bipartite and A(Gsb) < S, Gsb has an edge coloring 
that uses at most S colors. By Lemma [ 2 ] such a coloring defines a solution x' such that ([7]), ([8]> hold. ■ 

(Theorem [lj: We already showed that for bipartite networks JTK-MMK solves JTK and JTC-BIP solves 
JTC. Lemma [l] concludes the proof. ■ 
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(Lemma [5]l: Let z*,y* be an optimal solution for JTK. Without <[T|) and <[8j, JTIC would be equivalent to 
MMK with some restrictions on its parameters (unit-size MMK items and two MMK capacity values). Therefore, 
let z',y' be the solution returned by JTK-MMK, U(z',y') > aU(z*,y*). Finally, the solution is feasible due to 
Lemmas [3] and [4| ■ 

(Theorem [3): Let Bo be the set of BSs with no backhaul links (BSs whose degree is 0 in Gj). In line [2] of 
JTK-MAT the selected transmissions are determined using Ammk- This set of transmission is an a-approximation 
with respect to a JTK instance with B" = Bq. Therefore, to complete the proof we can assume that every BS has 
a backhaul link and show that the transmissions determined in line [7] are an (2a/3A(Gj))-approximation. 

Let Zq JS . 2 /) )ls be an optimal solution for OJS. The sum of weights for all edges in C, as computed in line [5] of 
JTK-MAT, is at least a£/(zo JS ,yojs)- 


Any graph G has an edge coloring using at most (|)A(G) colors 128 . Such a coloring for Gj partitions C 
into (|)A(Gj) matchings. Since in line [6] of JTK-MAT a maximum weight matching £ is obtained, the sum of 
weight for edges in £ is at least )a’(Gj) = (2a)/(3A(Gj))(7(Zojs,2/oj S ). Let z' ,y' be the solution returned 
by JTK-MAT. Then, U(z',y') > (2a)/(3A(Gj))U(z* oss ,y* OJS ). 

We now show that z',y r is feasible. Note that <[5j and ([6} hold since the solution associated with each edge {a, 6} 
is feasible (line |5| and solutions of different edges in £ use items whose h(i) is in disjoint (no two edges in £ 
share a vertex). Due to line [7] of Algorithm JTK-MAT, in the returned solution if z[ = 1 and h(i) = {a,b} then 
{a, b} £ £. Since £ is a matching, Gsb is bipartite and A(Gsb) < S- Using Lemma [4] we conclude that z' ,y’ is 
feasible. 

Finally, by Lemma [2] JTC-BIP solves JTC. By applying Lemma [I] we complete the proof. ■ 

(Theorem |4|: Let Zq [S be the optimal solution for OJS. The sum of weights for all vertices in B, as computed 
in line [5] of JTK-STA, is at least a:L r (z ( “ )ls ). In each iteration of the repeat loop (line [7] of JTK-STA), in line 10 


the 

// 


utility added to the solution z" equals U (6 max ) (note that z" gets updated to 1 at most once due to the update of T 
in line [12|. In line [12] 6 max and its neighbors are removed from consideration. The utility lost due to this removal 


is at most A(Gj)(7(6 max ). Therefore, the total utility of z" returned in line 18 is at least (a/A(G j))U(Zq ]S ). 

To show that z" returned by JTK-STA is feasible, it is sufficient to note that due to the correctness of A MMK 
function SOL-STAR returns a feasible instance with respect to I" and that I (line [ 9 ]) in different iterations contains 
packets of disjoint sets of transmitting BSs (therefore having positive weight only in disjoint sets of dimensions). 

Finally, by Lemma [2] JTC-BIP solves JTC. By applying Lemma [l] we complete the proof. ■ 

(Lemma |6j: It is clear that if Gj is planar, Gsb is also planar. Therefore, to complete the proof it is sufficient 
to show that Gsb is series-parallel. 

To show Gsb is series-parallel, we use the following definition | |35) . A multigraph is series-parallel if it has no 
subgraph isomorphic to a subdivision of a clique of size 4. Since Gj has no subgraph isomorphic to a subdivision 
of a clique of size 4, by adding parallel edges to Gj such a subgraph cannot be created in Gsb- Therefore, by the 
above definition Gsb is series-parallel. ■ 

(Theorem [2): The following result, mentioned in [36], is needed. 
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Let Eu C E denote edges in E whose both vertices are in U and let 5(G) = max{^j : U C V, 

\U\ > 3 and odd}. If G is planar and series-parallel then x'(G) = max{A(G), [5(G)]} and JTC-PSP 
from [361 finds an edge coloring that uses x'(G) colors. 

Recall that in JTK-PSP an instance for MMK with D+ B 0 a dimensions is constructed. The weight constraints 
for the new |£> 0 dd| dimensions are equivalent to requiring that for a feasible solution z', 5(Gsb) < S. Therefore, 
for such Gsb there exists an edge coloring that uses at most S colors, and by Lemma [2] such a coloring defines 
a solution x' such that |7]), ([8]> hold. Since JTK-PSP invokes Ammk which returns an a-approximation solution to 
the constructed MMK problem, JTK-PSP is an a-approximation for JTK. 

Finally, by Lemma [2] JTC-PSP solves JTC. By applying Lemma [T] we complete the proof. ■ 

(Theorem |sj: Let A € A. In order to demonstrate positive recurrence of {L(i)} t >o, we define a Lyapunov 
function, and show that it has negative drift outside some closed set of states. Let l = (Zi, h ,..., In, In) and define 
the quadratic Lyapunov function V(l) = Y^n =t In + l'n- 
Consider the one-slot drift 

AV(Z) = E {V(L(t + 1)) - V(L(t)) | L{t) = l}. (13) 

By Q and ((2} we compute 

L n (t + l) 2 = L n (t) 2 + 2 L n (t)(w n (t) - n^\t) - ^(t)) + (W n (t) - \t) - ^(t)) , 

L n (t + l) 2 = L n (t) 2 + 2L n (t)(^\t) - - 2) (t )) + (^ ] (t) - MnH*)) , 

Substituting this into ( f]~3j ) we obtain 

N 2 2 

AV(l)=J2n(w n (t)-n { n 1) (t)^^(t)) +($Ht)-$Ht)) \L{t)=l} 

n—1 

N 

+ 2Y j nLn(t)(w n (t)-ti\t)-^\t)) + L n (t)(^(t)-^\t)) | L(t) = I}. (14) 

n—1 

Since the W n (t) have finite second moment, and the fin\t) have finite support, we can bound (for some constant 

C < oo), 

N 2 2 

E { {Wnit) - ^\t) - ^\t )) + (vffit) - ^ 2) (i)) I L(t) =1} < C. 

n= 1 

The second part of ( fl4| ) can be written as 

N 

2 Y j nLu(t)(w n (t) - nPit) - ^(t)) + L n (t)(^\t) - fi) 2 \t )) I Lit) = 1} 

n—1 

N 

= 2 (inK - l n n^n\t) + ^(t) I Lit) =1}- i n E{^\t) - ^(t) I Lit) = /}) . (15) 

n—1 

Since A £ A, we know by the definition of A that there exists a flow vector / such that the conditions in © 
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hold. Moreover, A £ A also implies that there exists are Conv(i?) such that fn^ if fn ^ > 0, j = 1,2,3, 

n = 1,..., N. Since / is dominated by r, and r £ Conv(f?), there exist <ti, ..., ctr such that 


1*1 


1*1 



f — y&i r ii y ~ < i, 

*=i *=i 

(16) 

where r, = (rj^, rfj , ,.. 

Using (fl6| we obtain 

., r l'%) represents the i-th vector in R. 



1*1 

In An = Zn(/l 1} + f^) = L ^3 ( r i,n + r S >) ! 

2=1 

(17) 


1*1 

o = L(f^ 2) - /F) = L E ^ ( r ln - • 

(18) 


i=l 


By combining ( fT7) > and |I8}, and exploiting the structure of the fj, chosen according to uq (see ( fl2| )) 

N N |*| 

»EU-!EE« (Wrg + r«») + US - rS)) 

n—1 n= 1 2—1 

JV 1*1 

< 2^ (z n (E{ M W(f) + M ( 3 )(f) I L(t) =1}) +i„(E{^(t) -^ 3) W | £(i) = I})) £ a,. (19) 

n=l 2—1 

Substituting ( fl9l ) into (13) yields 

N 

2 ^ E{L„(i) (w„(i) - - /4 3) (*)) + L n {t)[^\t) - /# \t )) | L(i) = Z} 

n—1 


|-R| AT 


< - 2 1 - 


2(i - E a ‘) E (w^ 1} (*) + z4? } (i) I i(i) = 1} - Z n E{^ 3) (*) - ^u\t) I L(t) = Z}) < 0, (20) 

2—1 n=l 

where the last inquality follows from the choice of /.in 1 m The expression in ( [20] ) can be made arbitrarily small 
by increasing the state Z. Positive recurrence of {L(t)}t>o then follows from 114 Theorem 2.2.4], ■ 
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Appendix B 
Extended Results 


In this appendix we relax assumptions (i)-(iii) from Section III and show the applicability of our algorithms 
to the case where a packet can be transmitted using one of several Modulation and Coding Schemes (MCSs). 

we 


This appendix contains three sections. In Section B-A we present the extended network model. In Section B-B 


introduce the extended OJS Problem and in Section B-C we describe the decomposition framework and extend our 
algorithms. 


A. Network Model 

Each wireless transmission requires an MCS, which is selected from a finite set Ad = {1,..., M} of supported 
MCSs. We define the configuration set 1Z = {0} U M. such that r > 1 indicates wireless transmission using MCS 
r, and r = 0 indicates a packet forwarded to another BS over the backhaul. Thus, a pair (i,r) (% £ I. r £ 7Z) 
defines a wireless transmission of packet i with MCS r (r > 1) or packet i forwarded over the backhaul (r = 0). 
Note that the pair (z,0) is not feasible if f3(i) = 1. 

For each packet-configuration pair (i,r), we define two properties: its size T(i, r) and its success probability 
p(i,r). In case of forwarding over the backhaul, I’(«, 0) represents the size of packet i in bytes. For r > 1, the 
r(r, c) represents the number of scheduled blocks required for transmission with MCS r, which depends on the 
packet size in bytes T(i, 0) and the MCS r. 

The success probability p(i, r ) represents the probability that packet i will be successfully received by the user 
n(i) (r > 1) or that packet i forwarded over the backhaul is successfully received by BS(n) (r = 0). We assume 
that p(i,m) is higher if f3(i) = 1 compared to fi(i) = 0, since the SINR of a user is greater when a packet is 
joint-transmitted. The capacity of a backhaul link between BSs a and b is l a b bytes. 


B. The OFDMA Joint Scheduling (OJS) Problem 


We now formulate the extended version of the joint scheduling problem presented in Section IV Capacity 
constraints apply both to the subframes of the B BSs as well as the C backhaul links, and we denote the total 
number of such constraints by D = B + C. We order these constraints such that the constraint b corresponds to the 
BS b , b = 1..... /I, and constraints B + 1,..., D correspond to the backhaul. Define the D-dimensional capacity 
vector K = ( K\ ,..., Kd) such that for 1 < d < B, I\d = S (number of scheduled blocks), and l(,i = l a i, for 
d>B + 1. 

In order to describe the capacity used by the wireless transmission of packet i, we introduce 


h(i) 


{BS(n(*))| if P{i) = 0, 

{BS(n(*)),BS(?r(i))} if fi(i) = 1. 


If (3{i) = 0 then h returns only the serving BS, and if 3(i) = 1 it returns both the serving and secondary BS. 
We define the function w : T x 1Z —» (Nq) 15 , where w(i,r ) denotes the D-dimensional vector that represents the 
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capacity used for (i,r). If r > 1, then w(i,r) is the all-zero vector except for [w(i,r)]b = r (z, r) for b G h(i). 
If r = 0, then only the entry corresponding to the appropriate backhaul link is positive, and equal to the length 
in bytes of packet i. The function u : X / TZ K>• R + represents the utility of scheduling packet i according to 


configuration r, and is essentially an extension of the utility function presented in Section IV 


Based on the set of packets 1 and the utility function u, the centralized scheduler determines the set of 
transmissions to take place in the upcoming subframe, by maximizing the aggregate utility. The scheduler must 
also determine which scheduled blocks will be used for each packet transmission, such that for JT the scheduled 
blocks of the serving BS and secondary BS are aligned (i.e., have identical index). The scheduling decisions are 
represented using indicator variables z ir G {0,1} and x irs G {0,1}, where z ir indicates if a transmission (z,r) 
takes place, and Xi rs indicates if scheduled block s is used by a transmission (i, r). The scheduler needs to solve 
the following integer programming problem (with z = (z ir ) ie x,reTi and x = ( x irs)iex,rentes)- 
OFDMA Joint Scheduling (OJS) Problem: 

I R 

^2^2z ir u(i,r) =: U(z) 


max 

X'Z 


2—1 r= 1 


T—1 

VI 

> 

Viel, 

(21) 

r—0 

I R 

[i,r)\ d <K d , VI <d<D, 

(22) 

2=1 r=0 

S 

^ ^ %irs — %ir 

S = 1 

■T(z, r), Vz G I Vr > 1, 

(23) 

R 

E E 

Xi rs < 1, Vh G B Vs G S, 

(24) 

r= 1 { 2 : beh(i)} 

Zir ^ {0? 1}? 

ViGlVrGK, 

(25) 

%irs ^ 

MiGlMr gKMsGS. 

(26) 


Constraint ( |2Tj ) ensures a packet is scheduled only once; ( |22| i ensures capacities in each subframe and backhaul 
link are not exceeded; ( [23} ensures sufficient scheduled blocks are allocated for each (z, r) transmission; and ([24} 
ensures that each scheduled block is used at most once in the subframe of each BS. 


C. OJS Problem - Algorithms 


1) Decomposition Framework: Similar to what was done in Section IV for the basic problem, we can decompose 
OJS into two parts as follows. 

Joint Transmission Knapsack (JTK) Problem: 


max U(z) 

Z 

s.t. (|2T), 3x: ([23), ([24}, @ hold 




TABLE III 

The performance and input required for different algorithms within 

THE DECOMPOSITION FRAMEWORK [Aj TK , Aj TC ] 


TABLE IV 

Algorithms for MMK 
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^4-mmk 

Ratio 

Tmmk(I, R, D, S) 

DP (2\L 

Optimal 

0(S u IRD) 

PTAS 1261 

1/(1 + «) 

0((IRf D/c '<) 

Greedy ||21 

OO 

0(1 Rlog(I R)) 


^4 JTK 

^4 JTC 

Ratio 

Running time 

Input Gj 

JTK-MMK 

JTC-BIP 

ty 

0(T M mk(I,R,D,S)) 

bipartite 

JTK-MAT 

JTC-BIP 

3A (Gj) 

0(CT M uk(I.R,3,S)) 

any 

JTK-STA 

JTC-BIP 

*(Gj) 

0(B j T mmk (I, R, 2A(Gj) + 1, 5)) 

any 

JTK-PSP 

JTC-PSP 

Ol 

0(2" + T M uk(I,R,2 ±1 ,S)) 

planar ser.paral. 


Joint Transmission Coloring (JTC) Problem: 


given z, find x s.t. ([23]), <[24]), ([26] hold 


Solving JTK to obtain z and then solving JTC to obtain x, provides a solution to OJS. 

The scheduled block graph for this version of JTC is the same as what was presented in Definition [2] except 
that every pair (*,r) such that r > 1 and Z{ r = 1 contributes T(r, c) edges to Esb instead of 1. 

2) Algorithms for General Graphs: We now describe the algorithms JTK-MAT and JTK-STA for the extended 
model. The performance ratio of JTK-MAT and JTK-STA in the extended model are the same as those in Section [y] 
The proofs use similar ideas to those used in Section [V] and are omitted for brevity. Tables III and IV summarize the 
different options for solving OJS using the decomposition framework. Finally, note that the stability result presented 
in Theorem [5] can also be extended to this setting. 


Algorithm JTK-MAT Based on matching 
1: for b e B do 

2: If b has no backhaul link, run A MMK to solve a JTK instance with B' = {6}, I' = {i £ I : h(i) = B'}, 

C = {} and set z ir as determined by A MMK . 

3: for e = {a, 6} G C do 

4: Run A mmk to solve a JTK instance with B' = {a, b}, I' = {i £ 1: 3r € 1Z, h(i , r) C B'} 

5: Assign U(z) found by Ammk as a weight for e 

6: Compute maximum weight matching on Gj and store the result in the edge set 8 
1: Set Zi r = 1 only if Be G 8 such that z rr = 1 in the solution returned in line [4] for edge e 
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Algorithm JTK-STA Based on star subgraphs 
l: function SOL-STAR(6, W, C , T) 

2: Run Ammk to solve a JTK instance defined with B' = {b} U {a : {a, 6} G C'}, C = {{a, b} : {a, b} G 

C', a G B 1 , be B'}, T = {i G T : h(i) C B'} 

return z as determined by Ammk 

3: for h G B do 

4: Z[b\ g- SOL-STAR(X, B , C, X) 

5: Assign £/(.Z[&]) as a weight for b in B 

6: Initialize £" g- £>; C" g- C; X" g- X 

7: repeat 

8: Find the vertex fo max in £1" with maximum weight 

9: X G- {i G I" : 3r, 2[& m ax]t;r = 1} 

10: for all i G X, z" G- Z[6 max ] ir 

11: J G- {a G B" : 36 G B", {& max , 6} G C", {a, 6} G C "} 

12: X" G- X" \ X; Z3" G- B" \ {a : {a,6 max }GC"} 

13: Remove from C" edges with an endpoint not in B" 

14: for b G J do 

15: Z[b\ G- SOL-STAR(X, B", C", X") 

16: Update U(Z[b]) as a weight for b in B" 

17: until B" is empty 

18: return z" 





